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Abstract 

Received-energy test for non-coherent decision fusion over a Rayleigh fading muhiple access channel 
(MAC) without diversity was recently shown to be optimum in the case of conditionally mutually 
independent and identically distributed (i.i.d.) sensor decisions under specific conditions HI, ||2l. Here, 
we provide a twofold generalization, allowing sensors to be non identical on one hand and introducing 
diversity on the other hand. Along with the derivation, we provide also a general tool to verify optimality 
of the the received energy test in scenarios with correlated sensor decisions. Finally, we derive an 
analytical expression of the effect of the diversity on the large-system performances, under both individual 
and total power constraints. 

I. Introduction 

Starting from classical distributed detection [3l, large efforts in the recent literature have been devoted 
to the implementation of distributed detection in wireless sensor networks (WSNs) Q, Q, IH. Local 
decisions in a WSN are usually transmitted to a decision fusion center (DFC) in order to improve 
reliability of geographically distributed sensing through central processing. Common system architectures 
make reference to the availability of parallel (non-interfering) channels from the sensors to the DFC fTl, 
m, However, more sophisticated setups have been investigated, where the intrinsically interfering 
nature of the wireless channel is exploited and not combated HI, ifTOl . ifTTl . 

Recently, in IH and 121, the received-energy test was studied for non-coherent decision fusion over a 
multiple access channel (MAC). More specifically, in |1] the received energy was claimed as optimal for 
the no-diversity case with conditionally (given the phenomenon) mutually independent and identically 
distributed (i.i.d.) sensor local decisions, as long as the probability of false alarm of the generic sensor 
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is lower than the corresponding probabiUty of detection. Also, analytical performances of the received- 
energy test in the diversity scenario were derived. However, optimality property of the test was not 
investigated. The optimality of the test for the no-diversity case with conditionally i.i.d. sensor local 
decisions was proven in |l2l|. Only the case with sensors whose probability of false alarm is lower than 
the corresponding probability of detection was considered. Nonetheless, the diversity case was still ignored 
in the optimality analysis. 

The main contributions of this correspondence are: 

• a rigorous proof of the optimality of the received-energy test for non-coherent decision fusiorQ over a 
Rayleigh fading MAC with arbitrary order of diversity and with conditionally mutually independent 
but non identically distributed (i.n.i.d.) sensor decisions, as long as each sensor probability of false 
alarm is lower than the correspondent probability of detection; 

• as a side result, a sufficient condition on the log-likelihood ratio (LLR) of the number of active 
sensors suited for testing received energy optimality in scenarios with con^elated local decisions; 

• analytical derivation of large-system performances for conditionally i.i.d. sensor local decisions as 
a function of the order of diversity, where two different scenarios are considered: (a) sensors with 
an individual power constraint (IPC); (b) sensors with a total power constraint (TPC). 

It is worth noting that in [11] a different scenario was analyzed, where: (i) conditionally i.i.d. sensor 
decisions were considered, and (ii) instantaneous channel state information (CSI) at the DFC was 
assumed. The focus was on the performance analysis of several sub-optrmal fusion rules in terms of 
complexity, required knowledge, probability of detection and false alarm. 

The paper is organized as follows: Sec. HIl introduces the system model; in Sec.|lll]we present the main 
results of this correspondence, while in Sec. |IV] we draw some concluding remarks; proofs are confined 
to Appendices. 

Notation - Lower-case (resp. Upper-case) bold letters denote vectors (resp. matrices), with a„ (resp. 
fflra,™) representing the nth (resp. the (n, m)th) element of the vector a (resp. matrix A); upper-case 
calligraphic letters, e.g. A, denote discrete and finite sets; denotes the N x N identity matrix; On 
(resp. Iat) denotes the null (resp. ones) vector of length A^; E{-}, (•)*, (•)"'', K(-), and ||-|| denote 
expectation, transpose, conjugate transpose, real part, imaginary part and Frobenius norm operators; P(-) 

'Although, energy receiver and non-coherent are not synonyms, in the paper we will confuse them. In the related literature, 
such a misuse is common due to the fact that the energy detector is the default receiver adopted for non-coherent decision 
fusion. 
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and p{-) are used to denote probability mass functions (pmf) and probability density functions (pdf), 
while P{-\-) and p{-\-) their corresponding conditional counterparts; A/'c(/^, S) (resp. AA(/i, S)) denotes 
a circular symmetric complex (resp. real) normal distribution with mean vector and covariance matrix 
Xl, B{k,p) denotes a binomial distribution of k trials with probability of success p and Xl denotes a 
chi-square distribution with L degrees of freedom; {a*b){£) denotes the convolution between series a{i) 
and b{i); finally the symbols ~ and A- mean "distributed as" and "convergence in distribution". 

II. System Model 

A. WSN modeling 

We consider a distributed binary hypothesis test, where K sensors are used to discriminate between 
the hypotheses of the set T-L = {Ho,Hi}, representing, not necessarily, the absence (Hq) or the presence 
(Hi) of a specific target of interest. The A;th sensor, /c € /C = {1, 2, . . . , K}, takes a binary local decision 
dk € Ti about the observed phenomenon on the basis of its own measurements. 

Each decision d^ is mapped to a symbol G = {0, 1} representing an On-Off Keying (OOK) 
modulation: without loss of generality we assume that dk = Hi maps into Xk = i, i ^ {0, 1}. The quality 
of the kth sensor decisions is characterized by the conditional probabilities P{xk\Hj). More specifically, 
we denote Po.k — P {x^ = and Ppk = P {xk = ^Hq), respectively the probability of detection 

and false alarm of the A;th sensor. 

The sensors communicate with the DFC over a wireless fiat-fading MAC, modeled through i.i.d. 
Rayleigh fading coefficients with equal mean power. The DFC employs an A^-diversity approach in 
order to combat signal attenuation due to small-scale fading of the wireless medium. The diversity can 
be accomplished with time, frequency, code or antenna diversity (as recently proposed in |[TOl . |[T2|| ). 
Statistical CSI is assumed at the DFC, i.e. only the pdf of each fading coefficient is available. 

We denote: y„ the received signal at the nth diversity branch of the DFC after matched filtering and 
sampling; hn^k ~ M: (O, o"^) the fading coefficient between the fcth sensor and the nth diversity branch 
of the DFCo; Wn the additive white Gaussian noise at the nth diversity branch of the DFC. The vector 

^It is worth noting that assuming an asymmetric model for channel coefficient statistics would be more realistic. However, this 
would make the results much more dependent on the specific scenario without adding any significant insight from a theoretical 
point of view. A symmetric model for channel coefficient statistics is assumed for a two-fold reason: on one side it can be 
considered as a starting point before analyzing more realistic application-dependent scenarios; on the other side a symmetric 
scenario could represent scenarios in which power control is considered. 
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model at the DFC is the following: 

y = Hx + w (1) 

where y G C^, H G C^^^, x G Af^, ii; ~ A/'c(OAr, (t^/at) are the received signal vector, the channel 
matrix, the transmitted signal vector and the noise vector, respectively. Finally, we define the random 
variable i = i{x) = J2k=i ^k, representing the number of active sensors and the set £ = {0, ... , K} of 
possible realizations of £. 

B. LLR 

The optimal test 131, fl3\ for the considered problem can be formulated as 

H=Hi 



A. 



opt 



In 



' p{y\Hi 
.piy\Ho)\ 



7 



(2) 



H=Ho 



where H, Aopt and 7 denote the estimated hypothesis, the Log-Likelihood-Ratio (LLR, i.e. the optimal 
fusion rule) and the threshold to which the LLR is compared to. The threshold 7 can be determined 
to assure a fixed system false-alarm rate (Neyman-Pearson approach) or can be chosen to minimize the 
probability of error (Bayes approach) Q, lfT3l . An explicit expression of the LLR from Eq. Q is given 
by 



A, 



opt 



In 



in 



1 



\\y\\\^ 



PU\H, 



sr^ 1 



\yr- 



PU\Ha 



(3) 



where we have exploited the conditional independence of y from (given i). 

In the case of conditionally (given Hi) i.i.d. sensor decisions i{PD,k, PF,k) = {Pd,Pf), A; G /C) 
we have that i\Hi ~ B{K,P£)) and i\HQ ~ B{K,Pp). Differently, when local sensor decisions are 
conditionally i.n.i.d. the pmfs P{£\Hi) are represented by the more general Poisson-Binomial distribution 
lfT4l . ifTSl . |[T6l . with expressions given by: 

K K 

- \{(PD,kr\{{i-PD,s)^^~''-^ (4) 

x:l{x)=tk=l 



P{i\Hi) 

pmo) 



K 



s=l 
K 



(5) 



x:e{x)=ik=l s=l 

It is worth noting that Eq. ^ requires sums which are infeasible to compute in practice unless the 
number of sensors K is small. For this reason different methods have been proposed in literature for its 
efficient evaluation. The alternatives include fast convolution of individual Bernoulli pmfs |[T4l . recursive 
approaches |[T5l and a Discrete Fourier Transform (DFT) based computation 111 61 . 
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III. Optimality of Received Energy 

As already stated in [1], the received energy ip = ||y||^ is a sufficient statistic for the LLR, since 
Eq. Q depends on y only through i/j. However, sufficiency alone does not guarantee that the test 

H=Hi 

^ i (6) 

H=Hq 

is equivalent to Eq. As shown in 0, the test in Eq. Q is optimal iff Aoptitp) is a strictly increasing 
function of ip. If this property is satisfied, the test in Eq. ([2]) is equivalent to Eq. Q by simply setting 
7' = A~p-i-{'y). For this purpose in the following we first introduce a general optimality test (in the form 
of a sufficient condition) which relates the pmfs P(i\Hi), Hi G Ti, to assure that Aoptitp) is strictly 
increasing in the case of an A^-diversity MAC. 

Proposition 1. A sufficient condition for Aopt{ip) to be a strictly increasing ffinction of ^ is given by: 

A(^)>A(£-1), £eC\{0} (7) 

where X{£) = In 



P{t\Ho) 



Proof: The proof is reported in Appendix [A] ■ 
The above proposition states that strictly increasing property of assures optimality of the test in 
Eq. We will refer to A(^) as the LLR hereinaftefl It is worth noting that such a sufficient condition 
is independent of the order of diversity N. Also, Eq. ([7]) depends on the WSN model only through the 
number of active sensors I (given Hi) and does not require any specific assumption on P{x\Hi), e.g. 
conditional mutual independence of local sensor decisions, i.e. P{x\Hi) = Hfe^i P{xk\Hi). This means 
that Eq. (Q plays the role of a general property for received energy optimality, to be verified even in the 
case of conditionally correlated local sensor decisions. 

In the simplest case of conditionally i.i.d. local sensor decisions, (PD.k, Pp.k) = {Pd,Pf), G /C, as 
assumed in HI, 121, the strictly increasing property of A £, is equivalent to 

(f)P|,(l-P^)^-^ {e^)P'n\l-PD)''-'^' 



)P|(1 - Pf)^-^ {,-i)Pf\^ - Pf)^-'+' 

that reduces to Pd > Pp. This result, not only confirms theoretical findings for optimality of Eq. ^ 

when = 1 as in [Tl, 13, but it also proves the optimality of the test over the Diversity MAC (i.e. 

^Note that we will not consider in Eq. Q (and throughout the paper) the case I — Q when testing ^-LLR strictly increasing 
property, since A( — 1) has no physical meaning. 
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> 1) used in 111; this result shows the effectiveness and simplicity of Proposition [T] w.r.t. the approach 
taken in 13. 

Differently, when sensor decisions are conditionally i.n.i.d. (i.e. the case of a heterogeneous WSN), 
the following theorem generalizes the result in Eq. dU. 

Theorem 1. If P{x\Hi) = [IfcLi ^(^fcl^i)' G V., and Po^k > PF,k, k £ IC, the i-LLR satisfies the 
strictly increasing property described in Eq. ([7|) and thus Eq. ([6]) is the optimal test when an N -diversity 
MAC is employed. 

Proof: The proof is reported in Appendix |Bj ■ 
Theorem [T] states that under non identical sensors {PD,k, PF,k) and A^-diversity MAC the received 
energy ip is again the optimal test. Also, Theorem [J relies on a sufficient condition, i.e. specific WSN 
configurations not satisfying the assumption Po.fe > PF,k, k £ IC, but still verifying Eq. d?]) may exist. 
Although such a case is not "typically" of interest, since the condition PD,k < PF,k for kth sensor is not 
realistic in practical scenarios (i.e. sensors operating under nominal conditions), it proves the robustness 
of the received energy in scenarios with some faulty (or byzantine) sensor^. 

We finally evaluate analytically the performances as the number of sensors goes large, in the case 
of conditionally i.i.d. sensors. Both IPC and TPC on the WSN and arbitrary diversity N are considered 
here. This result generalizes [2], where no-diversity (N = 1) and IPC assumptions were made in deriving 
formulas. We define 

A 1 /"Hx 



^ = > ^ + ^ (9) 



where, compared to Eq. (O, 2 ^ merely scaling factor and Hx is replaced with ^= in z in 

order to keep a fixed amount of average energy e = E{\\Hx\\ } w.r.t. A^. Then we define the system 
probabilities of false alarm and detection, respectively Pp^ and Pp^^, as: 

PFo=P{\\zf>l\Ho), PDo=Pi\\zf>l\Hi), (10) 
Eqs. ^ and (fTOl) hold for TPC scenario when replacing z with z = . ^ ^ ( + ) . In this case 

\'Pfo-^ \VKN ) 

the average energy is kept fixed w.r.t. both K and N . 

"♦a WSN with K sensors such that [Pd.i,Pf,i) = (0.5,0.05), (Po,2,Pf,2) =(0.4,0.1) and (Pd,3,Pf,3) = (0.3,0.4) 
verifies the property in Eq. Q. 
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Theorem 2. If af^ and cx^ are finite, as K ^ +00 

z|Ho^Ac(oA,,-ijA,^, z\HiANt:{oN.-^lNy (11) 

w/iere = p^^i+a^ N '^"'^ = p^f2!|^^2 jy - ^^^'^ large-system {Pdo-ipc,^Fo-ipc) '^^^ S^ven by: 



Pk-iPcil) = exp(-7iV)x J]l(7Arr 



(13) 




Pk-lPc{l) = exp(-^) X >^3(^) ; (14) 
while {Ph^-TPcPpo-TPc) given by: 



AT-l 




(15) 



PDo-TPci^) = exp(-^^) X >^-(^^) . (16) 

Proof: The proof is reported in Appendix O for the IPC case; performance in the TPC scenario can 
be derived in a similar fashion. ■ 

As expected, if = 1 the result of Eqs. (fT3l) and ([141) coincides with the one given in (T. Sec. IV]. 

It is worth remarking that, in both IPC and TPC scenarios with diversity, Neyman-Pearson and Bayesian 
error exponents are zero (cfr. with |2]), because the large-system ROC can not be driven toward the 
point {P^^,PpJ = (1,0) by increasing the number of sensors, as long as the diversity A^ is kept 
finite. This intuition is confirmed by the non-zero values assumed under an IPC and a TPC by the 
large-system J-Divergence, J{p{z\Ho),p{z\Hi)) = N x + -21, J{p{z\Ho),p{z\Hi)) = 



N X 



yPpari 



I PpCXp 



, which represents a lower-bound for the system error probability \YT\, 



thus enforcing a zero Bayesian error exponent. Differently, the Neyman-Pearson error exponent is given 
by lim;^^+oo -i^^i^^^^^, under Pf.{i,K) < a. If we choose 7^ such that Pp^_^{'ya) = a, then 
limx->+oo - In [1 - Pooila, K)] = - In [l - P^oila)] < +00, giving again a zero error exponent. 

Note that the performance in TPC scenarios differ through the ratio {ap /ao) < 1 (cfr. Eqs. (fT3] ) and 
([141 ) with Eqs. (fTSl ) and ([T6l )) which represents the performance reduction factor w.r.t. IPC scenarios. 
Note that apjaiy: (i) is an increasing function of the ratio ^ (i.e. the received SNR), with limiting 
value equal to one; (ii) is a decreasing function of A^, meaning a diverging separation in performance 
between IPC and TPC as A^ increases. 
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The diversity affects in a significant way the large-system probabiUties of detection and false alarm, 
under an IPC, by shifting the Receiver Operating Characteristic (ROC) toward the upper-left comer, as 
shown in Fig. [T] meaning a performance improvement. Differently, it can be seen how a different effect is 
present in the TPC case, where an increase of N does not always coincide with performance improvement, 
but rather an optimal N, depending on {Pd^Pf, ^), exists. This effect was already noticed in HI and 
it is due to non-coherent combining loss of branch contributions. 

Finally, in Figs. |2] and [3] we verify, through simulations, the convergence of the ROC to the large system 
expression {K — > +00) given by Eqs. ([T3] ) and ([141 ) (resp. Eqs. ([TSl l and ([T6l)). under IPC (resp. TPC). 
It is apparent that the convergence under the TPC is faster w.r.t. the IPC case, because in both cases the 
large system ROC expressions rely on the Gaussian approximation of the Gaussian mixture given by Eq. 
([3]). For such a reason, for a given K, imposing a TPC on the WSN assures a better matching w.r.t. to 
the IPC case, since all the components of the mixture will be more concentrated. 

IV. Concluding Remarks 

In this correspondence we showed the optimality of the received-energy test for decision fusion 
performed over a non-coherent diversity MAC with conditionally i.n.i.d. sensor decisions. We derived 
a sufficient condition on the LLR of the number of active sensors which can be applied to test the 
received-energy optimality in WSN with conditionally correlated sensor decisions. Finally, we showed, 
through analytical results, how the diversity in a WSN with conditionally i.i.d sensor decisions affects 
the large-system performance under both IPC and TPC. 
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■□- ' ROC - K = 500 
□ ' ' ROC - K = too 
B - ROC - K - 50 



Figure 2. ROC comparison: large system vs finite number of sensors (K £ {50,100,500)) under IPC. £ {lj2}, 
(a"h/o"m)ds = 15. 




Figure 3. ROC comparison: large system vs finite number of sensors (K G {50,100,500)) under TPC. TV G {1,2}, 
icrl/al,)dB = 15. 
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Appendix A 
Proof of Proposition 1 

We prove in this section that Eq. O is a sufficient condition for the optimality of t/j test. From Eq. ([3]), 
looking at the LLR as a function of V', we get: 



dip 



1 



1 

a(V') 



■ K 

E 

ii=0 
K 



dtp 
dg{il;,£2 



K 



£2=0 



5-0 



K 

E 

£1=0 



(17) 



where we denoted g{Tp, i) 



'4> 



(-g.2 +£o-^)iv exp (^— ^2 ^f^i J and a{ip) indicates a positive function of ^ (i.e. 
a(^/^) > 0, VV' G M+) . Strictly increasing property of LLR is guai-anteed if ^^"q^"^^ > 0, V0 G M+, thus 



manipulations from Eq. (fTTl ) lead to 

£1-1 



i,/,^ ,^ ^ r^9('0,^i) / , . N dg{i;,£2) , 



Oil) 



> 



(18) 



£1=1^2=0 

where A;(^i,^2) = [P(^i|i^i)P(^2|^o) - P{i2\Hi)P{ii\Ho)]. In deriving Eq. ^ we could express the 
double sums in Eq. (fTTl) as a function only of the indices £1 > £2, since the term in bracket in Eq. ([TSl l 
equals to zero when £1 = £2- Noting that ^^gj^^'^ = — ^£f^2-^ g{ip, £) the condition is rewritten as 



K 



YYk{£i,£2)9{i^,£i)g{i^,£2) 



<{h-£2) 



i^l+hal){al+£2alj, 



> 



(19) 



Since both g{il^,£) and the term in square brackets are positive (note that indices in the sums are such 
that £1 > £2), the term k{£i,£2) is responsible for the sign of each term in the sum. Then a sufficient 
condition for Eq. ( fT9l ) is obtained assuming that each of those terms is positive. This is achieved if the 
following property holds 



ki£i 



> 0, 



h > £2. 



(20) 



It is easy to demonstrate that the condition k{£,£ — l) > 0, £ G £\{0}, representing the strictly increasing 



property of £-LLR, i.e. 



P(i\Ha) 



P{i-l\Ha) 



A^H , is equivalent to Eq. (I20l ). In fact Eq. (I20l ) implies that 



LLR is strictly increasing; this is verified just substituting 



?i — 1. Differently, we can show that 



^-LLR strictly increasing property implies Eq. (l20l ) by constructing the chain of inequalities p|^|j^^| 



> 



P(£i-l\Hi) 
P{Ii-1\Hq) 



> 



> 



P{i2\Hi) 

P(e2\Ho) 



■jM, all deriving from £-LLR strictly increasing property. 
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Appendix B 
Proof of Theorem [H 

We prove the strictly increasing property of ^-LLR by induction. Let us assume there exists a set of 
{t — 1) sensors with local performances {Po^k, Pp.k), A; G {1, . . . , t — 1}. The number of active sensors 
in this case is denoted it-i = X]fc=i ^fc' ^t-i G ^t-i- We denote the probability of t active sensors 
out-of-(t — \), given Hi, as i-t„i(£|i/j), Hi £ Ti and the corresponding £-LLR as Aj_i(^). 

Initialization: the strictly increasing property of ^-LLR in single sensor case Ai(^i) > Xi{ii — 1), 
ii € £i\{0}, is straightly verified when Po i > Pf,i- 

Induction: Let us assume that for a specific configuration of {t — I) sensors the ^-LLR Xt-i{it-i) 
satisfies the strictly increasing property, that is Xt-i{it-i) > Xt-i{it-i — 1), it-i G ^t-i\{0}. If we add 
the tth sensor satisfying Pd ^ > Pp^t and we prove that the new £-LLR Xt{it) > Xt{it — 1)> G ^t\{0}, 
i.e. it retains strictly increasing property, then the proof is complete. 

To proceed let us first define a{i) = Pt^i{£\Hi), b{i) = Pt-i{i\Ho), c{£) = Pi{£\Hi) and d{l) = 

Pi{e\Ho). 

The number of sensors transmitting when the tth sensor is added is then given by it = X]fc=i = 
It-i + xt. The pmfs Pt{(t\Ho) and Pt{£t\Hi) are then given by HH 

PmHo) = {b * d){£t) Pt{it\Hi) = {a * c){£t) (21) 
The LLR strictly increasing condition is then formulated as follows 

{h*d)[£t) {b*d){£t-l) 
By exploiting the support set of c{£) and d{£) we can rewrite Eq. (l22l ) as follows 

Efcg{o,i}c(fe)«(^t -fe) Efcg{o,i}c(fc)a(^f-l-fc) 
Zke{o,i} mm -k)^ Efc6{o,i} d{kWt -l-k) 

where obviously a{t) = b{t) = 0. Exploiting c(0) + c(l) = (1 - Pd,*) + Pd,* = 1, d{0) + (i(l) = 

(1 — PF,t) + = 1, we obtain 

[l-c{lM£t) + c{l)a{£t-l) [l-c{lM£t-l) + c{l)a{£t-2) 
[1 - d{l)]b{£t) + d{l)b{£t - I) [1 - d{l)]b{£t - 1) + d{l)b{£t - 2) 

The condition expressed in Eq. (I24l ) can be rewritten as: 



{[1 - c(l)][l - d{l)][a{£t)bi£t - 1) - aiit - l)6(^t)]}+{c(l)d(l)[a(4 - l)b{£t - 2) - a(£t - 2)6(4 - 1)]} + 

+{c(l)[l - d{l)][ai£t - l)b{£t - 1) - aiit - 2)6(4)]}-{[l - c(l)](i(l)[a(£t - l)b{£t - 1) - a(4)6(£t - 2)]} > 

(25) 
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Since fgj > > we have that: 

[aiet)b{£t - 1) - a(^t - > [a{£t - - 2) - a{it - 2)b{£t - 1)] > (26) 

[a(^OK^t - 2) - a(^t - 2)b{it)] > (27) 

The condition in Eq. ( [25] ) is satisfied since positivity of the first two terms follows from the inequalities in 
Eq. (l26l) . and the difference of the third and fourth terms in Eq. (1251 ) is positive because c(l)[l — d{l)] > 
[1 — c(l)]d(l) (since Pd ^ > PF,t) and exploiting the inequality in Eq. (l27l) . This concludes the proof. 



Appendix C 
Proof of Theorem [2] 

The proof follows in the first part similar steps as in 121; for this reason we will only sketch it 
and underline the substantial differences. We use here the characteristic function of the vector z\Hi, 
i e {1,2}, denoted as <I>2|//, to easily evaluate the limit for K — )• +oo. We then use this result, 
in conjunction with Levy's Continuity Theorem |18] to demonstrate the convergence in distribution of 



large-system p{z\Hi). Let us now write the characteristic function of z\Ho as a function of t = {t\,t^2 

^z\Ho{t) = ^z\Ho{eMjt{zi + jtiz2)} = /•••/ exp{jt\zi + jtiz2) X ^p{zi, Z2\e)P{e\Ho)dzidZ2 

(28) 

where zi = 'R{z} and Z2 = (with tj, i G {1,2}, representing the index-corresponding dual 

vectors over Fourier domain). Following analogous steps as in [2|, exploiting: i) conditional independence 
assumptions such as p{zi,Z2\t) = p{zi\£)p{z2\i) and p{zi\£) = Y[^=iP{^i,s\^)y * £ {1,2}; ii) the 
characteristic function of x ~ M{0, cr^) is given by <I>x(t) = exp(- 



[181; we get 



K 



1 ^ 



exp 



s=l 



1 ^ 

s=l 



+ 



cr:, 



NKPf (TiKPF 



Pp exp 



1 ^ 



2,s) 



s=l 



NKPf 



K 



+ {1-Pf) 



(29) 



where in the last line we exploited £\Hq ^ B{K, Pp). Also, exploiting similar noteworthy limits as in 



[[3, eventually we have that limx- 



exp 



1 v^AT 

2 Z^s=l 



2N 



Applying the Continuity 



Theorem flS], we obtain z\Hq — )• A/c (On, j^In) ■ In a similar way it can be shown that z\Hi 
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The last part consists in proving Eqs. ([T3] ) and ([T4l) . The large-system probabilities of false alarm and 
detection can be expressed in the equivalent form: 

P*fM = mi^f >7l^o) = p(^C>7l^o) (30) 

P*dM = P{\\zf>l\Hi) = p(^-^p^C>j\Hi^ (31) 

where ^ ~ X^(2N)- ^^'^ probabilities are then easily calculated evaluating the cumulative distribution 
function of ^ [1811: 




which provides the result. 
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